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We introduce a general formalism to analyze nucleation phenomena in inhomogeneous media which 
considers the influence of the metastable phase, which is treated as a heat bath in which clusters 
are embedded, in the dynamics of the nucleation process. A kinetic equation for the evolution of 
the probability density of the clusters, taking into account possible inhomogeneities in the bath, 
is derived using mesoscopic nonequilibrium thermodynamics. From it, we obtain the nucleation 
rate and analyze the role played by the metastable phase in the process. We discuss in particular 
condensation and polymer crystallization in a temperature gradient. 

I, INTRODUCTION 

Nucleation is the initial step in many phase transformations, and thus plays a crucial role in many processes of 
scientific and technological importanceii2iiLiiti. Despite many years of intense investigation, it is still a very challenging 
topic due to the inherent difficulties associated with the nonequilibrium nature of the phenomenon, and to the fact 
that the most interesting entity, namely, the embryo of the new phase or nucleus, is a unstable and highly improbable 
object. Due to these difficulties, the analysis has been mostly focused on the simplest case of homogeneous, isothermal 
and isotropic nucleation. But the real situation under which nucleation very often occurs does not correspond to those 
conditions. One has to take into account that the real process occurs in a medium, which in general has spatial, thermal 
or velocity inhomogeneities, and which in turn may exert a significant influence in the process. 

One of the situations in which this fact becomes more evident is in polymer crystallization. Contrary to the case of 
solidification of simple substances, the crystallization of polymers occurs over a wide range of temperatures. Moreover, 
in realistic situations, the external conditions, especially those concerning the temperature, may change with time and 
thus alter the process significantly. These characteristics make it necessary to study the process under nonisothermal 
and inhomogeneous conditions. 

Our aim in this paper is to develop a more realistic model of nucleation and crystallization that considers the 
possible influence of the medium. We will focus on the simplest cases in which the medium may affect the kinetics. 
These cases involve the presence of spatial inhomogeneities in the thermodynamic properties or in the conditions 
under which nucleation occurs; the presence of temperature gradients; and, finally, the influence of flow and stress in 
the nucleation process, a feature that will be discussed in the accompanying paper (henceforth referred to as Paper 
II). 

In the case of homogeneous isotropic nucleation, we can leave spatial dependencies aside since the process occurs 
identically at any point within the system. We can thus focus our description on the evolution of cluster size as 
a function of external conditions (pressure, temperature, density). Homogeneity thus justifies the use of a global 
thermodynamic approach. 

However, when the system is inhomogeneous, the conditions controlling condensation and crystallization vary from 
point to point and therefore a local description of the process must be considered. 

For this purpose, we will divide the whole sample into volume elements small enough so that, within them, temper- 
ature, pressure, density, etc. may be considered homogeneous, but large enough to guarantee that thermodynamics 
can still be applied. In terms of length scales, this means that the typical length of the volume element must be much 
smaller than the gradients present at the system and much larger than any length at the molecular scale. 

It is important to remark, that nucleation involves two clearly different length scales of respective macroscopic 
and mesoscopic natures. The macroscopic scale is the one on which thermodynamic quantities such as pressure, 
temperature, density, etc. vary, and these quantities can be considered as uniform within each elementary cell, even 
though they may change along the sample. On the other hand, nucleation occurs on a mesoscopic scale. 

For that reason, mesoscopic nonequilibrium thermodynamics (MNET) 6 - 1 * becomes the natural framework within 
which analyze the process. MNET is a powerful, and systematic method to describe the kinetics of nonequilibrium 
processes occurring on the mesoscopic scale. In its present formulationi^, it constitutes a generalization of the standard 
nonequilibrium thermodynamics, capable of describing the kinetics of mesoscopic entities in terms of arbitrary variables 
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of state or degrees of freedom. It has been applied successfully to many out-of-equilibrium processes, and also to the 
kinetics of homogeneous nucleationiiifiiii. 

The key ingredient in extending the formalism to the present situation (i.e. an inhomogeneous medium) is to 
carry out a local description resorting to a reference state in which local equilibrium^ is assumed to hold. Starting 
from this reference state we use the MNET formalism in order to obtain the Fokker-Planck equation describing the 
dynamics of the nucleation process in this inhomogeneous medium. The resulting Fokker-Planck equation provides a 
complete description of the system at the mesoscopic level. However, it contains more information than is accessible 
from experiments, performed on a macroscopic scale. We then describe the macroscopic evolution of the system using 
hydrodynamic equations, constructed by averaging the relevant quantities with the probability density occurring in 
the Fokker-Planck equation. 

In this way, one obtains the kinetic equations governing both nucleation and transport in the system. The formalism 
is able to recover some results of the kinetic theory, but in a more simple and systematic way, within a nonequilibrium 
thermodynamics framework. In particular, we obtain relaxation equations that can describe the short-time inertial 
regime of the dynamics. However, for nucleation, this short-time regime is usually not relevant, hence we focus our 
analysis on a long-time regime where diffusion is the dominant mechanism of transport. 

For simplicity, we will not consider angular variables. These can be easily incorporated into our scheme, but their 
inclusion complicates the description of the process. For the sake of concreteness, we investigate in detail two inter- 
esting cases. In the first situation, analyzed in the present paper, we derive the kinetic equations governing nucleation 
in spatially inhomogeneous systems under the presence of a temperature gradient. The second situation, discussed 
separately in Paper II, involves nucleation in a stationary flow, representative of the influence of stresses or mechanical 
processing. These examples resemble situations of practical interest. Using the results obtained within our theoretical 
framework, we finally discuss their implications for both condensation experiments and polymer crystallization^. 

II. NUCLEATION UNDER SPATIALLY INHOMOGENEOUS NON-ISOTHERMAL CONDITIONS 

Consider a nucleation process taking place in a non-isothermal situation where the metastable phase, globally at 
rest, is subjected to conditions that create a stationary temperature profile T(x), where x is the position coordinate. 

The presence of spatial inhomogeneities makes it necessary to perform a local description in terms of /(n, x, u, t) /N, 
the probability density of finding a cluster of n G (n,n + dn) molecules at position x e (x,x + dx), with velocity 
u € (u, u + du), at time t. N is the total number of clusters in the system. The whole system is then characterized 
thermodynamically by the local energy density per cluster e(x,i), the entropy density per cluster s(x,i), and the 
number density 



where p m is the density of the heat bath, assumed constant. 

Entropy variations of the system, can be expressed in terms of the corresponding local thermodynamic quantities 
through the Gibbs equation (see Appendix of Paper II) 



where p(n, x, u, t) is a generalized chemical potential. Notice that the last term in the previous equation is reminiscent 
of the corresponding one for a mixture in which the different components would be specified by the pair of indices 
n and u. The generalized chemical potential can be specified through the use of the Gibbs' entropy postulate. This 
formula establishes the connection between statistical mechanics and thermodynamics through the definition of the 
entropy in terms of a probability distribution function, i.e. 



where fi eq and psi eq are, respectively, the probability distribution and the entropy at local equilibrium (indicated by 
the subscript leq). Variation of this expression yields 




(1) 




(2) 




(3) 




(4) 
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where the variation of the local equilibrium entropy satisfies the fundamental equation of thermodynamics 



5{ps)leq = y<5(pe)/ eg - pieqSpieq, (5) 

and (j,i eq (x) denotes the local equilibrium chemical potential which is independent of n and u , but in general will 
depend on position. By comparison of the previous equation with the Gibbs equation one obtains the expression 
for the generalized chemical potential 



p(n, u,x,t) = k B T\n- h/Xi eg (x), (6) 

Jleq 

essentially the chemical potential for a species identified by the indices n and u. 

At local equilibrium, the system, viewed as a gas of clusters in the metastable phase, is described at each position 
by an equilibrium distribution function of the internal variables (i.e velocities and sizes) in which the temperature 
is that at x. Consequently, the local equilibrium distribution is just the probability of occurrence of an equilibrium 
fluctuation represented by a cluster of n particles and velocity uSfl 



t i i\ ( Pleg - C(n,U) \ 

Jleq{n,u,x,t) = exp I ^-^ j. (7) 

The quantity C(n, u) denotes the free energy cost of formation of a cluster of size n, velocity u and mass m(n) = m\n, 
and is given by 



C(n, u) = AG{n) + ^m(n)u 2 , (8) 

where AG(n) represents the free energy of formation of a cluster of size n at rest and the second term is its kinetic 
energy. By using the local equilibrium distribution 0, the chemical potential can be written as 

p,(n, x, u, t) = fc^Tln f(n, x, u, t) + C(n, u) , (9) 

which corresponds to that of an ideal system since interactions between clusters have been neglected. 

The next step is to formulate the conservation laws governing the evolution of the relevant quantities of the system, 
namely the probability density /, the internal energy e(x, t), and the entropy s(x,£). In absence of external forces, 
the continuity equation for / may in general be written as 

d f d t dJn nn\ 

W = - u ' v/ -^- J "-^T' (10) 

where J n is the flux of clusters in size space, and 3 U is a new current resulting from the interaction of the clusters 
with the heat bath. 

The conservation law for the local energy density is£*i& 



dpe 
~~dt 



= -V-J g , (11) 



and states that, in the absence of external forces and viscous heating effects, the total internal energy of the fluid 
element at x can only be altered by a heat flux J g ,. 

The rate of entropy change per unit volume can be evaluated by differentiating the Gibbs equation l|2} with respect 
to time. Thus 



= t^-7t 7^ l u-J-dndu. (12) 

dt T dt T J dt y ' 



Introducing the energy balance (1 1 H and the continuity equation IjlOJl . we obtain 

dt ~ T V Jq 



~ /(u-V/ + !rJ« + ^Wndu. (13, 



T I V du dn 
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Integrating by parts, assuming that the fluxes vanish at the boundaries, and introducing Eq. © for the chemical 
potential leads to the following balance equation for the entropy 



where the entropy flux 3 S is given by 



J s = ^ ^ - k B J u/(ln/ - l)dndu J fu {AG{n) T AH{n)) dndu , (15) 
and the entropy production a, which must be positive semidefinite according to the second law of thermodynamics, is 

<J = ~— j'-VT- — / J„- —dndu- — / J n — dndu. (16) 



T 2 q T J du T J n 3n 

In Eq. Illot . AH(n) is the enthalpy of an n— cluster and J' q is the irreversible heat flux 

J' q = 3 q - J fu ^AH{n) + im(n)u 2 ^ dndu , (17) 

in which the integral represents the transport of heat due to the diffusion of clusters. In order to obtain the entropy 
flux we employed the identity 



J V • (u/) In fdndu = V ■ (J u/(ln/ - 1) dnduj , 
and the thermodynamic Gibbs-Helmholtz relation 18 



(18) 



{AG(n)\ AH(n) . , 

V(^^MJ- ^VT, (19) 

valid when the system is in mechanical equilibrium. 

The entropy production can be interpreted as a sum of products of conjugate forces and fluxed. In the present 
situation, Eq. allow us to identify the thermodynamic forces 

1 „ m 1 du 1 du . . 

conjugate to the heat flux J 9 , velocity flux J M , and size flux J n , respectively. Thus it becomes clear that the origin of 
the thermodynamic forces are the presence of gradients or, in accordance with Eq. ©, variations of the distribution 
function with respect to its local equilibrium value. 

Now, according to the tenets of nonequilibrium thermodynamics, we may assume linear relationships between forces 
and fluxes. Assuming isotropy, which implies that quantities of different tensorial rank are not coupled, and locality 
in the internal space for which currents at each point (n, u) are determined only by the local properties, one obtains 

4 = -^f VT - J i L Tu ■ ^ dndu , (21) 
J^-^VT-Iwg, (22) 
/ -_ir ?E 

where by postulate the phenomenological coefficients satisfy the Onsager reciprocal relations, e.g. 

L T u = -LuT ■ (24) 
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It is useful to rewrite the phenomenological coefficients in a more convenient way. Defining A 



conductivity, D n 



— knL n 



f 



as the diffusion coefficient in n-space, and j3 



kff- as the thermal 
n jT" and £, = ^fjr as friction coefficients, 



and using the expression for the chemical potential, Eq. Q, the currents can be written as follows 



J' q = -AVT 



£m fu 



k B Tdf 
m du 



dndi 



(25) 



1 \ m du 



(26) 



/ --D \°L 

on 



1 dC 



(27) 



k B T dn ' 

Introducing these expressions into the continuity equation, Eq. ifTfil) . we finally obtain the Fokker-Planck equation 



on 



I dC 
k B T dn 



d_ 
du 



/3 /u + 



k B T8f 
m du 



(28) 



that governs the evolution of the inhomogeneous density distribution of clusters in a bath having a non-uniform 
temperature distribution. The transport coefficients (3 and £ still remain unspecified. Their identification must be 
carried out through a proper interpretation of the macroscopic relaxation equations derived from the Fokker-Planck 
equation. We discuss this issue in the next section. 



III. HOMOGENEOUS NUCLEATION IN SPATIALLY INHOMOGENEOUS SYSTEMS: DIFFUSION 

REGIME 



The Fokker-Planck equation, Eq. I|28J1 . provides the complete description of the evolution of the probability dis- 
tribution of clusters at mesoscopic level. In particular, it can even describe inertial regimes at the earlier stages of 
the nucleation process and retains information about the velocity of clusters. However, the process of equilibration in 
velocity space is usually much faster than those in size or coordinate space. After times much larger than the charac- 
teristic time for velocity relaxation, the system enters the diffusion and thermal diffusion regime where the evolution 
is governed by a simpler set of equations. In that regime, we can describe the system in terms of the moments of the 
distribution function, which are related to the hydrodynamic fields, namely the reduced probability density, defined 
as 



fc(n, x, t) = J f(n, x, u, t)du , (29) 

the velocity in (n, x)-space 



v c (n,x,t) = / ( T i / ufdu, (30) 

and the second moment 



V = J (u-v c )(u-v c )/du, (31) 



which corresponds to the kinetic part of the pressure tenso^iS* 



The evolution equations for these moments follow by taking time derivative of Eqs. ll29|l - ll3"T)l and using the Fokker- 
Planck equation <|2Bj! 9 " 1 ™- The result is the continuity equation 
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dt dn 
the equation for the conservation of momentum 



l> fc = - V-/cV c - ^- I J n du, (32) 



fj£ = -V • 7> - I3f c v c - |/ C VT - f(u - v c )^-J„du, (33) 
and the equation for the evolution of the pressure tensor 



^-V = -V • Q - 2(V ■ Vv c ) s -PV-Vo- - ^^/ c /3 
at m 



I 



d 

(u- v c )(u- v c ) — J„du. (34) 



In the previous expressions, Q = J(u — v c )(u — v c )(u — v c )/dv is related with the kinetic part of the heat flux, a 
superscript s means symmetric part of a tensor, and the total or mobile derivative is defined as 

In a similar way, we could derive the evolution equations for the higher-order moments of the distribution, which 
constitute a coupled hierarchy of hydrodynamic equations&i^Sii. The closure of the hierarchy can be performed 
through time-scale considerations^. 

The momentum conservation law 1(331) enables us to identify the forces acting on the cluster, and the transport 
coefficients f3 and £. The third term on the right of that equation represents a force acting on the cluster, arising 
from the presence of the temperature gradient, and giving rise to the thermophoretic effect or thermophoresis 24 - 2 ^. 
Its origin can be naively understood on microscopic basis^. The cluster immersed in a thermal gradient is being hit 
by particles of the host fluid. Since collisions with particles from the "hot" side are stronger than those from the "cold" 
side, a net force driving the particle toward the cold side appears. The force per unit mass on a particle is given by 

F T = -— VT, (36) 
m 

where ? is called the thermophoretic coefficient. By comparison with the third term on the right of Eq. (|33l) we can 
establish the following identification 



£ = ?T. (37) 

The second term on the right of Eq. ll33)l can be identified with the hydrodynamic force exerted by the fluid on the 
cluster, with j3 playing the role of the friction constant. This quantity establishes the characteristic relaxation time 
scale for the velocity. It can be estimated using Stokes' law, f3 ~ ^7^, which indicates that its magnitude is very 
large for small clusters. Consequently, the discussion of the behavior of the system may be carried out by expanding 
the hierarchy of evolution equations for the moments in powers of (3~ x . This time scale also motivates the separation 
of the dynamics into two well-differentiated regimes: an inertial regime for t -C characterized by the relaxation 
of the variables toward the diffusion regime, which is achieved when t ^> /3 _1 . 

A remarkable feature of our theory is the fact that it provides a simple and complete framework for accurately 
analyzing the dynamics of mesoscopic systems, even in inertial regimes at very short times. The set of equations 
113211 - l|34H which governs the hydrodynamic behavior of the 'gas' of clusters, are analogous to that obtained in kinetic 
theory2i)22. Notice that the continuity equation, Eq. IjlOJI . has a form similar to the Boltzmann transport equa- 
tion, with the flux J„ playing the role of the collision integral. Therefore, within the framework of nonequilibrium 
thermodynamics, our mesoscopic theory is able to reproduce the results of kinetic theory in a simpler way. 

For instance, by retaining subsequent orders in the expansion of the continuity equation, one can succes- 
sively recover corrections to the diffusion behavior, analogous to the Chapman-Enskog or Burnett corrections of the 
Boltzmann equationfi*^. 

The short-time inertial regime is usually not observed in typical experiments, focused on nucleation phenomena 
occurring on longer time scales. It is then sufficient to center the description in the diffusion regime and to truncate 
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the hierarchy of moments at the level of the pressure tensor, thus discarding contributions from the heat flux and 
higher order moments. We shall retain only the lowest order correction in /3 _1 . 

First the evolution equations for the pressure tensor, Eq. ll3"4T) . and the velocity, Eq. will be discussed within 

the diffusion regime. The terms appearing in the equation for the pressure tensor ll-'Ut . involve different time scales. 
In the diffusion regime, i.e. for t ^> P -1 , time derivatives can be neglected when compared with terms proportional 
to (3. Notice that, in this sense, the divergence term V • v c is essentially a time derivative (a fact that follows from 
mass conservation equation), and accordingly can be neglected. Moreover, velocity relaxation is usually faster than 
relaxation in the size-space. Consequently, we can also neglect the contribution arising from the current J n . Taking 
all these considerations into account, the equation for the pressure tensor then reduces to 

V = ^f c t, (38) 

TO 

which is that for an ideal gasfi^S, being 1 the unit tensor. Introducing the value of the pressure tensor into Eq. l(33jl 
for the evolution of the velocity yields 

fc ^t = ~ V ' (A)/3/cl) - ^^-/cVT - /3/ c v c - J (u - v c ) J- J„du, (39) 

where Dq — is the spatial diffusion coefficient. In the diffusion regime we can again safely neglect the time 
derivative and the term proportional to J n . The resulting equation can be written as 

Jd = / c v c = -V (A>/c) - ^VT/ C , (40) 

which defines the current in the diffusion regime. Inserting previous expression into the continuity equation ll32l) . 
yields 



^=V-(V( J D /c)+^- 1 VlnT/ c )-^ j J„du, (41) 
where the velocity integral of the current in n-space is 



We can reasonably assume that, D n , the diffusion coefficient in n-space is approximately independent of the velocity, 
but that in general the drift term does depend on u. In fact, from the energy barrier formula, Eq. @, this drift is 



8C 1 2 dAG(n) 

^ = 2 miU+ ^^- (43) 
However, in the diffusion regime the system will have achieved equilibration in velocity space. Therefore the relation 



^TOU 2 /du = ^k B Tf c 



(44) 



reminiscent of the energy equipartition, holds here as well as in equilibrium. The integral of the current in size-space 
can then be written as 



df c 1 ( dAG 3 k B T 
Tfo + ~k~B~f V dn + 2~^~ 

- D \dfc + J_dAG(j^ 
dn kgT dn 

where AG(n) is the height of a modified nucleation barrier whose value is 



(45) 
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AG(n) = AG(n) + -k B T\nn, 



(46) 



which includes in an averaged way the effects of the Brownian diffusion of clusters. This correction was discussed in 
more detail in Ref^, where rotational degrees of freedom were also taken into account. 

The evolution equation for the cluster distribution function in spatially inhomogeneous systems can then finally be 

expressed as 



df d 
-g = V • (V (Do/e) + f/T'V lnT/ c ) + — 



It is convenient to rewrite the spatial flux of clusters as 



1 dAG 



dn ksT dn 



fr 



(47) 



\ D = -D Vf c - ^-f c VT - ^—f c VT = -D Vf c - D th ^rf c , 
rap 1 1 



(48) 



where it is evident that Jd has two contributions: ordinary diffusion described by Fick's law and a drift term that 
can be identified with thermal diffusion. The quantity 



D th = D 



= A) 1 



is the thermal diffusion coefficient. Another interesting quantity is the thermal diffusion ratio 



(49) 



kr = (50) 

which measures the importance of thermal diffusion, also known as Soret effect. Finally, it is important to emphasize 
that the relation found in Refill between the thermophoresis and the Soret effect, 

? = M*r-l), ( 51 ) 
is automatically recovered, a fact that follows in a straightforward manner from Eqs. lIHTIl. (l49l and l(50j) . 



A. Temperature Evolution 

The above formalism not only prescribes the evolution of the probability density, but it can also describe the 
associated evolution of the temperature field that originates from the balance of internal energy. 

The local internal energy has two contributions: the internal energy of the clusters and that of the heat bath. If we 
neglect thermal expansion, variations of the total internal energy can be related to temperature variations through 
the thermodynamic relation pSe — cypST, where cy is the specific heat of the system, bath plus clusters, at constant 
volume. Neglecting viscous heating effects, the total internal energy may vary due to the presence of the heat flux 
and also due to the release of latent heat associated with the phase transition. The heat released in the formation of 
a cluster of the new stable phase of size n is £m(n)f c , where £ is the latent heat per unit mass. Consequently, the 
equation governing the evolution of the temperature field is 



-V§— + (52) 

where p f (x,i) = J m(n)f c dn represents the total density of the new phase (liquid in the case of condensation or 
crystallized material in the case of crystal nucleation) . 

We have seen that the coupling between thermal and diffusion effects modifies both the diffusion current and the 
heat flux. To show this feature it is more convenient to work with the unmodified heat current because it is usually 
the quantity measured in experiments 
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3 q = 3' q + J u{AH(n) + ^mu 2 )fdndu. (53) 



Substituting in this expression the equation for the flux, Eq. Ij25l) . and using 



2 

which is valid in the diffusion regime, we obtain 



VL—mu 2 fdndw. = I k B T3 odn, (54) 



J g = -AVT + J (AH(n) + k B T + £m)3 D dn. (55) 
Now, employing Eq. I|48J1 for the diffusion current 3 B , we finally arrive at 

3 q = — AVT - J D th k B T-pdn + J AH{n)3 D dn, (56) 



where 



is the modified heat conductivity. 

The equation for the evolution of the temperature can finally be expressed as 

cyp^ = V • (AVT) - V ■ J AH(n)3 D dn + V ■ J D th k B T^j^dn + (58) 

From it, we can identify the basic mechanisms responsible for temperature variations, namely heat conduction, 
convection, thermal diffusion, and the release of latent heat in the crystallization process, corresponding to the 
different terms on the right hand side of that equation, respectively. Notice that, by neglecting the contributions 
arising from the diffusion of clusters and the modification of the heat conductivity, one can recover the usual equation 
for the evolution of temperature field, namely 

Once having developed the general formalism describing nucleation in inhomogeneous media, our aim will be to 
analyze the influence that the presence of these inhomogeneities may play in the nucleation process. This is the topic 
of the next section. 



B. Influence of Diffusion and Thermal Diffusion in Nucleation Experiments 

1. Condensation Experiments in Thermal Diffusion Cloud Chambers 

Many systems in which nucleation occurs are spatially nonuniform. In fact, some experimental setups impose 
and take advantage of these inhomogeneities in order to bring about and measure nucleation rates. For instance, 
in nucleation experiments in diffusion chambers a temperature (thermal diffusion chambers) or a velocity gradient 
(laminar flow diffusion chambers) is imposed to generate a very narrow region in which supersaturation exceeds the 
critical value so that nucleation occurs. 

Different mechanisms exist through which the presence of inhomogeneities may influence the nucleation process. On 
one hand, gradients change the heat and mass transport in the system as well as the transport coefficients. Moreover, 
thermal and diffusion effects may induce convection, which strongly alters transport in the chamber. In that case, 
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the conditions of supersaturation and temperature under which nucleation occurs may not be accurately described 
in these experiments. This is an important problem that has been abundantly treated in the literatur e)^ and will 
not be discussed here. 

On the other hand, there is a direct effect due to the loss of sub-critical clusters through diffusion and thermal 
diffusion. If the region in which nucleation occurs is very narrow, sub-critical clusters may escape from that region 
before having had time to grow beyond the critical size. This mechanism may prevent nucleation or significantly 
reduce its rate. We will focus our analysis precisely on this effect, partially extending previous works by Becker and 
Reiss^i, and Shi et al*±. 

In real experiments in thermal diffusion cloud chambers, the imposition of a temperature gradient generates tem- 
perature, pressure and supersaturation profiles inside the chamber (for a representative illustration see Fig. 1 in 
Refi2&). Measurements are performed under steady state conditions, thus implying that both these profiles and the 
nucleation rate are stationary. A salient feature of this profiles is that, typically, the region in which nucleation takes 
place is quite narrow. The activated nature of nucleation implies that when supersaturation is below the critical value, 
nucleation is for all intents and purposes absent. 

Under these conditions, Eq. lUTI) can be simplified by assuming that nucleation takes place only in a small region 
of thickness do (see Fig. QJ. 

The loss of sub-critical clusters per unit of volume in that region by diffusion and thermal diffusion can be approx- 
imated by 



loss = -LJ V .J D d^^- = ±- {-DoVf c - D th Yff c ) . (60) 

The saturation ratio decays rapidly on both sides, and due to the fact that the concentration of clusters depends 
strongly on supersaturation, there are virtually no clusters outside that region. If the concentration of clusters is zero 
on the boundaries of the chamber, then one can estimate the gradients for diffusion as f c /d\ and (see Fig. 

for definitions of d\ and da) • Taking these considerations into account the rate of loss of n-sized clusters by diffusion 
and by thermal diffusion is 



- -JH D ofc - AhV mT-J-/ c , (61) 
dod do 

where d^ 1 = d^ 1 + ■ 

It is important to remark that these approximations for the diffusion and thermal diffusion terms of Eq. lUTjl are 
lower bounds that underestimate the actual value, since a sub-critical cluster will decompose immediately outside the 
region do, and thus the real gradients are consequently stronger. Introducing these simplifications and the stationary 
condition into Eq. l(T7jl . we obtain 

d ~ 1 „ 1 dhxT „ , x 

0=- F Jn-A)vr fc-D th - —f c , 62 

on ddo do dz 

where J n is the flux of clusters in size space. 

Since only order of magnitude estimates of diffusion and thermal diffusion effects are of interest, we adopt some 
additional simplifications. For simplicity, we will not include the nonequilibrium correction to the nucleation barrier 
given by Eq. Il4(it . This certainly modifies the actual value of the nucleation rate, but has a weak influence on the 
diffusion due to its logarithmic dependence on size. We will then adopt the CNT model for the flux of clusters^, 
namely 



an V/eg 

where f eq (z, t) is the equilibrium distribution of clusters given by 



J n = -k+(n)f eq ^- (A ) , (63) 



f eq (z,t) = N^expinlnS -8n 2/3 ), (64) 

where Ni is the total number of monomers; k + (n) is the rate at which an n-cluster gains monomers, 9 — -j^p is the 
non-dimensional surface tension, being a the bulk surface tension and Si the area per molecule of the bulk liquid; and 
S is the supersaturation. 
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It is also convenient to construct a dimensionless version of equation Ki2ll by resorting to time-scale considerations. 
It is evident that the effect of diffusion will be important when the time that the cluster takes to diffuse out of 
the nucleation region is smaller than the time it takes to grow beyond the critical size. The spatial and thermal 
diffusion coefficients, and the rate of growth of a cluster provide the proper time scales for each of these processes. 
The characteristic nucleation relaxation time (see Refp2&) is 

(65) 



2k+{n*Y 



where 5 = 3(rt*) 2 / 3 6 |_1 / 2 defines the width of the interval around the critical size n* = ( 3 f„ s ) 3 within which the 
height of the nucleation barrier has dropped by less than lfc^T. The diffusion and thermal diffusion time scales are 
in turn t b = n d f° \ , andm = dln y° , respectively. Therefore, the dimensionless parameters which control the 

^QK' 1 l dz ±J th 

influence of diffusion and thermal diffusion are a = — , and b = — . It is then clear, that a value of a or b larger 

tb ' Tth ' ° 

than unity implies that the time required for a cluster to nucleate is larger than the time it takes to move out of the 
nucleating region c?o by diffusion or thermal diffusion, respectively. In that case a reduction in the nucleation rate by 
loss of sub-critical clusters is expected. In contrast, values of a and b smaller than unity indicate that a cluster in the 
nucleation region has enough time to become stable by growing beyond the critical size before it diffuses away. In 
that situation, the nucleation rate is not significantly altered by diffusion or thermal diffusion. 

Introducing these parameters into Eq. I|62ll . and defining the reduced distribution function y = f c / f eg , and size 
x = n/n*, we obtain the dimensionless expression 



where e = S/n* . This equation is useful for estimating the effects of the inhomogeneities in the nucleation rate for 
different values of the parameters a and b. 

2. Numerical Estimates for the Case of Condensation 

To analyze the importance of diffusion and thermal diffusion in real nucleation experiments, equation 16611 has been 
solved numerically 37 - for different values of the parameters a and b. We have employed the usual boundary conditions 
y = 1 for n = 1, and y = for n — > oo. 

For the case of condensation, the expression for the diffusion and thermal diffusion coefficients can be borrowed 
from kinetic theory. In fact, the Brownian diffusion coefficient of an n— sized cluster is approximately given by2i 

D o( n > = — — w^r — ' ( 67 ) 

8P tot v27rM c si 

where Ptot is the total pressure and M c is the molecular mass of the host gas. The thermal diffusion coefficient can 
be approximated^! by 

D th = k T D(n) = a T xx(l - xi)D{n), (68) 

where or is the thermal diffusion factor, and x\ — p/Ptot is the mole fraction of the nucleating substance. One then 
obtains 

4 Peq SP tot s^M, 

and 



27tt (k B T) 



3 



Ap eq SP t ots\, 



=axxi(l — x\)ij), (70) 
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where (j) = (ddo)" 1 and V = iL jfd^ 1 . In the previous expressions, the parameters that can be changed significantly 
are the total pressure, and especially the equilibrium (coexistence) pressure which decreases exponentially with the 
inverse temperature. One can then infer that pressure is the main quantity that determines the importance of diffusion 
and thermal diffusion. At low enough pressures, a and 6 can be of the order of unity , thus indicating the potential 
relevance of diffusion and thermal diffusion under rarefied conditions. 

To be more precise, we discuss two particular situations. The first deals with the influence of diffusion and thermal 
diffusion in a real experiment 35 in a thermal diffusion cloud chamber, where n-pentanol is the nucleating substance. 
The second deals with the same problem using a model compound that imitates the properties of a typical sulfuric 
acid, a substance of great interest in atmospheric processes. Thermophysical properties of both substances are listed 
in Table ITvl Typical values of the geometric parameters in a diffusion cloud chamber are2&: d = O.lh, d% = 0.7h, and 
d 2 = 0.2/i, where h = 42.3 mm is the height of the chamber, which results in <f> = 3.6, and -0 ~ 0.1 for a representative 
value of 50K for the temperature difference between the upper and lower plates. 

Figures [2] and present the ratio between the rate of nucleation in the presence of diffusion and thermal diffusion 
effects and in their absence, as a function of the dimensionless parameters a and 6. Plots have been obtained from 
numerical resolution of Eq. (jfifil) using Mathematica, and the properties of n-pentanol at T = 260K , 5 = 9, and 
P tot = 30 kPaM.. 

The estimated values of the dimensionless parameters a and 6 corresponding to those experimental conditions are 
a = 5 10~ 8 and 6 = 8 10~ 12 . As one can observe from Figs. ElandEjl these small values of a and b imply that neither 
diffusion not thermal diffusion have a significant influence on the results of these particular experiments. This is 
mainly due to the high values of both the equilibrium and the total pressures under which these experiments are 
performed. 

The situation changes drastically when either the equilibrium pressure of the substance or the total pressure is low. 
To illustrate that fact, we have repeated the previous calculation using the model compound whose properties are 
listed in Tab. HH Fi gures and illustrate the influence of diffusion and thermal diffusion on the nucleation rate 
for this model compound at T = 293. 15K, S = 12. 5^ 4 i^ . When nucleation takes place at a normal pressure, e.g. 
Ptot = 1 atm, the resulting value for the parameter controlling diffusion effects is a = 1O~ 4 0, indicating again that 
diffusion is not especially significant. However, this model compound is peculiar in having a small saturation pressure. 
This indicates that if the total pressure is low enough, diffusion effects can be important. In fact, in experiments 
performed in thermal diffusion chambers, for reasons of stability the chamber must be operated at reduced pressure 
to avoid convectio ni? 1 "'? 2 . Using a realistic bound of P to t < 100p eg , one obtains that a > 3 for P tot < 10~ 2 Torr, thus 
implying a reduction of the nucleation rate by three orders of magnitude, as indicated in Fig. [2 This confirms that 
at reduced pressures diffusion may have a significant influence on the nucleation rate of substances. 

The properties of some pollutants commonly found at the atmosphere are similar to those of the model compound 3 ^. 
Therefore, the results of this section provide some insight into the possible effect that diffusion may have in in low 
pressure upper atmospheric processes. 

3. Condensation Experiments in Laminar Flow Diffusion Cloud Chambers 

The previous analysis remains valid for experiments performed in laminar flow diffusion cloud chamber a 39 41 . In 
these experiments, the imposition of a steady laminar flow is responsible for the appearance of a region in which the 
temperature abruptly drops and nucleation occurs. Once again, the zone in which nucleation takes place is narrow, 
and temperature gradients are also present. (Representative temperature, supersaturation, and nucleation rate profiles 
occurring in these devices are shown for instance in Fig. 6 of Ref4£). Therefore, we can use the results of the previous 
section to assess the significance of diffusion and thermal diffusion in these experiments. 

To be more precise, consider a typical experiment using n-pentanol as a nucleating substance. From Ref4i, repre- 
sentative conditions for these experiments are T = 260K, S = 10, and Ptot — 100 kPa. Introducing the properties 
of n-pentanol tabulated in Tab. lHvj) into Eqs. ll69ll and JZQl, yields the values a = 4 1O~ 9 0, and 6 = 8 10 -12 ^ 
for the dimensionless parameters. These very small values of a and 6 indicate that, under those conditions, these 
experiments in laminar flow diffusion cloud chambers are not significantly affected by the loss of subcritical clusters 
through diffusion or thermal diffusion. 

4- Influence of diffusion and thermal diffusion in Polymer Crystallization 

A simplified model of polymer crystallization can be proposed in terms of the coordinates of the center of mass 
of a spherical polymer crystallite and the number of monomers it contains. Under those assumptions, the analysis 
developed in the previous section remains valid for the description of the polymer crystallization process. 
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Although the theoretical description of condensation and crystallization can be similar, the crystallization process 
has some features that makes its description more difficult. For the case of crystallization, and especially for polymer 
crystallization, the identification of the transport coefficients is not so straightforward. Moreover, the attachment 
rate fc + (n) and especially the surface tension are not really known. For this reason it is very difficult to arrive at 
a detailed quantitative estimate of the diffusion and thermal diffusion effects in the case of polymer crystallization. 
Nevertheless, we present here a qualitative discussion of the probable influence of such effects. 

In the case of crystallization, the rate k + (n) at which units attach to the growing crystal embryo, is roughly 
proportional to the rate of diffusion of the units. This proportionality implies that any change in the rate of diffusion 
in the system alters the crystallization rate accordingly, since the mechanism of diffusion is also responsible for the 
growth of the crystals. The importance of this becomes evident for the case of crystallization under nonisothermal 
conditions. When a temperature gradient is present in the system, the diffusivity of monomers is changed by the 
appearance of thermal diffusion, as indicated by Eq. I|48JI . Moreover, the total diffusivity is no longer isotropic but is 
different in the direction determined by the temperature gradient. Since any reduction or increase of the diffusivity 
leads to the corresponding reduction /increase of the nucleation rate, this implies that the presence of temperature 
gradients may significantly alter the nucleation rate. We can estimate the potential magnitude of this effect. 

Thermal diffusion is proportional to the temperature gradient and to the thermal diffusion coefficient D t h, as 
indicated by Eq. 14811 . D t h is in general smaller than the normal diffusion coefficient in gases, but may become 
extremely significant in polymers. A measure of the importance of thermal diffusion follows from the ratio of thermal 
and diffusion coefficients, or Soret coefficient 



D Txi{l-xi)' y ' 

Typical values of S for gaseous and liquid mixtures range from 10~ 3 to 10~ 5 K~ x in order of magnitude, the smallness 
of which indicates that mass diffusion is dominant. For polymer solutions, however, the situation is completely 
different. The Soret coefficient increases with the molecular weight to the extent of becoming significant at high 
values of the molecular weight, as reported in Ref4^. Moreover, the small thermal conductivity of the melted polymer 
allows the release of latent heat involved in the crystallization process to set up very big temperature gradients, a 
result that follows from Eq. . Both factors imply that the magnitude of the thermal diffusion in polymers under 
non-isothermal conditions is significant, and therefore can influence the crystallization process importantly. 



IV. CONCLUSIONS 



In this paper we have outlined some of the principles of mesoscopic nonequilibrium thermodynamics (MNET) and 
have used this discipline to establish the basis for a complete description of the kinetics of nucleation in which the 
metastable phase acts as a heat bath in which the embryos of the new phase are embedded. We have explicitly shown 
that the presence of the heat bath may play a significant role in the nucleation process. 

To illustrate this, we discussed, in particular, the influence of inhomogeneities in the medium on the process of 
nucleation by analyzing in detail the process of nucleation in the presence of a temperature gradient. We derived 
a kinetic equation, of the Fokker-Planck type, that governed the evolution of the cluster size distribution, as it was 
coupled to the evolution of the heat bath. This equation was then used to evaluate the influence of the media in two 
particular situations of practical interest, namely, condensation nucleation experiments in diffusion chambers, and 
polymer crystallization. 

The results of our analysis suggest that, when experiments are performed at normal conditions, condensation in 
thermal and laminar flow diffusion cloud chambers is not significantly affected by diffusion and thermal diffusion 
effects. However, in rarefied media, as in upper atmosphere or for substances with low equilibrium vapor pressures, 
those mechanisms can become quite significant. 

Another situation in which the nucleation process can be drastically influenced by the presence of a temperature 
gradient occurs in the case of polymer crystallization. Polymers usually have low thermal conductivities, which favors 
the development of very large temperature gradients, and high values of the Soret coefficient. Since both factors control 
the importance of thermal diffusion, this is a clear sign that these effects may become crucial in polymer crystallization. 
Finally, it is worth mentioning that diffusion and thermal diffusion effects have proved to be important also in chemical 
vapor deposition, as analyzed in RefA 

The framework proposed in this paper can be extended to more general situations, such as the one in which the 
system is subjected to a shear flow. This situation will be treated in the accompanying paper, Paper II of this series. 
The method we have introduced thus goes beyond the classical formalism and can account systematically for the real 
conditions under which nucleation takes place. 
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Table I 



n-pentanol 



M = 88.15 g/mol 

p eq = 133.322 exp(90.079043 - 9788.384/T - 9.9 log T) Pa 
pi = 0.270 + 1.930229 Z 1 ' 3 - 8.414762 Z 2/3 + 19.226001 Z - 18.559303 Z 4/3 

+6.555718 Z 5/3 g/cm 3 
a = 26.85469 - 0.07889 (T - 273.15) dyn/cm 
M c = 4.0026 g/mol 

1/etT = (-0.7272 -T/(16.36- 0.2882 T))( Xl + 0.12281) +0.089303 



Model compound 




M -- 


= 94.13 g/mol 




Pi = 


-- 0.868 g/cm 3 (vi = 1.8 10" 3 cm 3 ) 




a = 


30.0 dyn/cm 




Peg 


= 10' 5 Torr 
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